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We investigate the relativistic equation of state of hadronic matter and quark-gluon plasma at 
finite temperature and baryon density in the framework of the nonextensive statistical mechanics, 
characterized by power-law quantum distributions. We study the phase transition from hadronic 
matter to quark-gluon plasma by requiring the Gibbs conditions on the global conservation of baryon 
number and electric charge fraction. We show that nonextensive statistical effects play a crucial 
role in the equation of state and in the formation of mixed phase also for small deviations from the 
standard Boltzmann- Gibbs statistics. 

I. INTRODUCTION 
(N 

^V The physics of high energy heavy ion collisions is a goldmine of problems in statistical mechanics and thermody- 

namics due to a large average number of particles involved and possible phase transition phenomena in the hot and 
dense fireball created during the collisions [1]. In relativistic heavy ion collisions the baryon density can reach values 
of a few times the saturation nuclear density and/or high temperatures. Furthermore, the future CBM (Compressed 
Baryonic Matter) experiment of FAIR (Facility of Antiproton and Ion Research) at GSI Darmstadt, will make possible 
to create compressed baryonic matter with a high net baryon density 2-4]. In this direction interesting results have 

i— i. been obtained at low SPS energy and are foreseen at a low-energy scan at RHIC [5r|9J. 

r~ Lattice calculations predict a critical phase transition from hadronic matter to quark-gluon plasma (QGP) at 

temperature T c of about 170 MeV, corresponding to a critical energy density e c w 1 GeV/fm 3 [T(|. On the other 
hand, in dense nuclear matter, baryons are forced to stay so close one to another that they would overlap. At 
large densities, constituent quarks are shared by neighboring baryons and should eventually become mobile over a 
C . distance larger than the typical size of one baryon. This means that quarks become deconfined and that at large 
densities and/or high temperatures they are the real degrees of freedom of strongly interacting matter instead of 

^-H ■ baryons. The process of deconfinement and the EOS of quark-gluon matter can in principle be described by quantum 

^ ' chromodynamics. However, in energy density range reached in relativistic heavy-ion collisions, non-perturbative effects 

in the complex theory of QCD are not negligible [lfj. The generated QGP in the early stages of the collisions does not 

^T ■ at all resemble a quasi- ideal gas of quarks and gluons because strongly dynamical correlations are present, including 

long-range interactions [l|, [ll|-fl4|. In the absence of a converging method to approach QCD at finite density one 

. ■ often turns to (effective) model investigations. Various results from QCD inspired models indicate that, increasing 

the baryon chemical potential in the phase diagram, a region of non-singular but rapid cross-over of thermodynamic 

observable around a quasi-critical temperature, leads to a critical endpoint (CEP), beyond which the system shows a 

first order phase transition from confined to deconfined matter. The existence or exclusion of a CEP has not yet been 

confirmed by QCD lattice simulations. Actually, there are some extrapolation techniques to finite chemical potentials 

15], although the precise location of the CEP is still a matter of debate [la |. For example, in Ref. [17] . the authors 

estimate the values T CEF = 162 MeV and // CEP = 360 MeV. Such a CEP can be in principle detected in future 

high-energy compressed nuclear matter experiments. 

Recently, there is an increasing evidence that the nonextensive statistical mechanics, proposed by Tsallis, can 
be considered as an appropriate basis to deal with physical phenomena where strong dynamical correlations, long- 
range interactions and microscopic memory effects take place [l8l - [2l| . A considerable variety of physical applications 
involve a quantitative agreement between experimental data and theoretical models based on Tsallis thermostatistics. 
In particular, in the last years there is a growing interest to high energy physics applications of nonextensive statistics 
|22j and several authors have outlined the possibility that experimental observations in relativistic heavy ion collisions 
can reflect nonextensive statistical behaviors [23TJ35J ] . 

The existence of nonextensive statistical effects should strongly affects the finite temperature and density nuclear 
equation of state (EOS) [36-38J. In fact, by varying temperature and density, the EOS reflects in terms of the 
macroscopic thermodynamical variables the microscopic interactions of the different phases of nuclear matter. The 
extraction of information about the EOS at different densities and temperatures by means of heavy ion collisions is a 
very difficult task and can be realized only indirectly by comparing the experimental data with different theoretical 
models, such as, for example, fluid-dynamical models [39). Related to this aspect, it is relevant to observe that a 
relativistic kinetic nonextensive theory [40| and a nonextensive version of a hydrodynamic model for multiparticle 
production processes have been proposed [41(. Very recently, nonextensive statistical effects on the hadronic EOS 
have been investigated by means of a Walecka type relativistic mean field model J42j. Furthermore, a nonextensive 
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version of Nambu-Jona-Lasinio model [43[ and the effects on color superconducting phase for two quark flavors due 
to a change to Tsallis statistics have been studied [44] , 

The main goal of this paper is to study how nonextensive statistical effects influence, from a phenomenological 
point of view, the nuclear EOS and, as a consequence, the relative phase transition at finite temperature and density 
reachable in high-energy heavy-ion collisions. Focusing our investigation to lower temperatures and higher baryon 
chemical potentials than the corresponding CEP values, a mixed phase of hadrons, quarks and gluons can be formed 
following the Gibbs conditions for the phase equilibrium. If, in general, we consider a substance composed of two 
conserved " charges" , like the baryon number and the isospin charge in heavy ion collisions at finite baryon density, 
the ratio between the two charges is fixed only as long as the system remains in one of the two pure phases. In the 
mixed phase, the concentration in each of the regions of one phase or the other may be different. Their values are 
restricted only by the conservation on the total charge numbers. The essential point is that conservation laws in 
chemical thermodynamics are global, not local. The main result of this formalism is that, different from the so-called 
Maxwell construction, the pressure in the mixed phase is not constant and therefore the nuclear incomprcssibility, for 
example, does not vanish [45j |. 

Furthermore, the scenario we are going to explore corresponds to the situation realized in heavy ion collisions 
experiments at not too high energy where finite temperature and high compressed baryon density is reached. In 
this condition, a not large fraction of strangeness can be produced and, therefore, we will limit ourselves to study 
the deconfinement transition from hadronic matter into up and down quark matter [46l448l | . We aspect that, in the 
range of temperature and density considered, the presence of strange particles does not significantly affect the main 
conclusions regarding the relevance of nonextensive statistical effects on the nuclear EOS. 

The paper is organized as follows. In Section [Til we introduce the basic formalism of the nonextensive statistics. In 
Section IlIIl we study the nonextensive hadronic EOS for symmetric and asymmetric nuclear matter and we explore 
the behavior of meson fields in presence of small deviations from the standard statistics. In Section ITvT we investigate 
nonextensive proprieties of the quark-gluon EOS. In Section [Vj we study the hadron to quark-gluon phase transition 
and the consequent formation of a mixed phase, mainly focusing our study in the variation of the first critical transition 
density for various set of parameters. Finally, we summarize our conclusions in Section [VII 



II. BASIC ASSUMPTIONS IN NONEXTENSIVE STATISTICS 

Nonextensive statistical mechanics introduced by Tsallis is a generalization of the common Boltzmann- Gibbs sta- 
tistical mechanics jla420| . It is based upon the introduction of the following entropy 

s q [f] = ^- l (i-J[f(x)]un\, (7/(x)dfi = iV (i) 

where /(x) stands for a normalized probability distribution, x and dQ denoting, respectively, a generic point and the 
volume element in the corresponding phase space. Here and in the following we set the Boltzmann and the Planck 
constant equal to unity. The real parameter q determines the degree of non-additivity exhibited by the entropy form 

©• 

The generalized entropy has the usual properties of positivity, cquiprobability, concavity and irreversibility, preserves 
the whole mathematical structure of thermodynamics (Legendre transformations). In the limit q —¥ 1, the entropic 
form (fTJ becomes additive and reduces to the standard Boltzmann-Gibbs entropy 

Si=- /"/(x)ln/(x)ffi. (2) 

Peculiarity of the Tsallis generalized thermostatistics is that if we have two statistically independent subsystems 
A e B, described, respectively, by the individual probability density f^ and f^ and we call f( A+B \x.A,*-B) = 
f( A ) (xa) f^ (x_b) the joint probability density of a composite system A+B, the nonadditivc (nonextensive) character 
of S q is summarized in the relation [20| 

S q [f {A+B) ] = S g [/ (A) ] + S q [fW]) + (I - q)S q [fM] S q [f^} . (3) 

In the limit q — > 1 , the third term in right hand side of Eq.dHJ) vanishes and the above equation reduces to the well- 
known additivity (extensivity) relation of the Boltzmann-Gibbs logarithmic entropy. Here, the word nonextensive 
should be associated with the fact that the total energy of long-range-interacting mechanical systems is nonextensive, 
in contrast with the case of short-range-interacting systems, whose total energy is extensive in the thermodynamical 
sense J20|. 



Second crucial assumption on nonextensive statistics is the introduction of the g-mean value (or escort mean value) 
of a physical observable A(x) 



, M(x) [/(x)]«d n 
/[/(x)]«dn 



w« = J T;, r^ • ( 4 ) 



The probability distribution can be obtained maximizing the measure S q under appropriate constraints related 
to the previous definition of the g-mean value. In this context, it is important to observe that the Tsallis classical 
distribution can be seen as a superposition of Boltzmann distributions with different temperatures which have a mean 
value corresponding to the temperature appearing in the Tsallis distribution. The nonextensive q parameter is related 
to the fluctuation in the temperature and describes the spread around the average value of the Boltzmann temperature 

Following the above prescriptions, it is possible to obtain the associate quantum mean occupation number of 
particles species i in a grand canonical ensemble. For a dilute gas of particles and for small deviations from the 
standard statistics (q « 1,) it can be written as [49l.[5C| 
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where (3 — 1/T and the sign (+1) is for fermions and (—1) for bosons. Naturally, for q — > 1 the above quantum 
distribution reduces to the standard Fermi-Dirac and Bose-Einstein distribution. Let us observe that nonextensive 
statistical effects vanishes approaching to zero temperature. This is the reason for which nonextensive effects could 
be significantly relevant in high energy heavy ion collision and probably, in the protoneutron star. In addition, in 
high density quark-gluon matter the color magnetic field remains unscreened (in leading order) and long-range color 
magnetic interaction should be present at any finite temperature, thus QGP appears to be an ideal candidate for 
finding some nonextensive behavior. 

Finally, let us observe that when the entropic q parameter is smaller than one, the above distribution have a natural 
high energy cut-off which implies that the energy tail is depleted; when q is greater than one, the cut-off is absent and 
the energy tail of the particle distribution (for fermions and bosons) is enhanced. Hence the nonextensive statistics 
entails a sensible difference of the power-law particle distribution shape in the high energy region with respect to 
the standard statistics. In this context, it is relevant to observe that in Ref. j43|, the authors postulate a modified 
quantum distribution function for fermions and bosons at the scope of satisfy the particle-hole symmetry, both for 
q > I and q < 1. In the present work we will focus our study for small deviations from the standard statistics and 
for values of q > 1, because these values were obtained in several phenomenological studies in high energy heavy ion 
collisions (see, for example, Ref.s [2J, [33 - t35| ). We have explicitly verified that, for the values q > 1 considered in 
this investigation, the prescription introduced in Ref. |43| does not affect the results that we are going to obtain and, 
therefore, we adopt the original formulation of nonextensive statistics. Furthermore, it is proper to remember that in 
a relativistic mean field theory, considered in this investigation, baryons are assumed as Dirac quasiparticles moving 
in classical meson fields, the field operator are replaced by their expectation values and the contributions coming from 
the Dirac sea are neglected. 

III. NONEXTENSIVE HADRONIC EQUATION OF STATE 

In this Section we study the nonextensive hadronic EOS in the framework of a relativistic mean field theory in 
which nucleons interact through the nuclear force mediated by the exchange of virtual isoscalar-scalar (a) , isoscalar- 
vector (to) and isovector- vector (p) meson fields |5ll453J |. As quoted in the Introduction, a similar approach has been 
studied in Ref. [421 ] for pure neutron matter and for symmetric nuclear matter (thus, without considering the effects 
of the p meson field) focusing principally the attention to a different range of density and temperature considered in 
this paper. Here, we are going to study the hadronic EOS at the scope of emphasize several features previously not 
investigated that result to be crucial for our following purposes. 

The Lagrangian density describing hadronic matter can be written as 

£ = £-QHD + £qfm , (6) 



where [53 



Cqhd = #7^ - (M - g a a) - ft^w" - g p ^f ■ p^ + -(d^ad^a - mla 2 ) 



\mlu^ + l -m 2 p p^ ■ p» - \f mv F>» - it 



-U(a) + -mlw^w" + -m 2 p p^ ■ ?» - 7 V" - jG^G^ , (7) 



and M = 939 MeV is the vacuum baryon mass. The field strength tensors for the vector mesons are given by the 
usual expressions F pv = dpLo v — d^uj^, G^ u = d p p v — d v p^ and U(a) is a nonlinear potential of a meson 

U(a) = ^aa 3 + ^ba 4 , (8) 

usually introduced to achieve a reasonable compression modulus for equilibrium nuclear matter. 

Following Ref.s [5J, [55[ , £ q f m in Eq. © is related to a (quasi) free gas of pions with an effective chemical potential 
(see below for details). 

The field equations in a mean field approximation are 

(t 7M a" - (M - g a a) - g^u - g pl °r 3 p)^ = , (9) 

m 2 a a + aa 2 + bcr 3 = g a < tptjj >= g G p s , (10) 

mluj = g u < -07°?/; >= g u p B , (11) 

™ 2 pP = 9 P < ^7°r 3 ip >= g P Pi , (12) 

where a = (a), u> = (io a ) and p = (p®) are the nonvanishing expectation values of meson fields, pi is the total isospin 
density, ps and ps are the baryon density and the baryon scalar density, respectively. They are given by 

P B = 2 E J 7^l n *(k)-Mk)}, (13) 

/A3h. a/t* 
p^WW+#)], (14) 

Z — Ti^p l 

where rii(k) and rii(k) are the q-deformed fermion particle and antiparticle distributions given in Eq.([5]); more explic- 
itly, in this context, we have 

ni{k) = [l + (q-l)p{E*(k)-p*)Yl^)+l ' (15) 

Mk) = ll + (q-l)P(E*(k)+ti)]i/(f-i)+l ■ (16) 



The nucleon effective energy is defined as Ei*(k) = y k 2 + Mi* , where Mi* = Mi — g a a. The effective chemical 
potentials p* are given in terms of the meson fields as follows 

(4 = IM, - 9u>u - T 3 tg P p , (17) 

where pi are the thermodynamical chemical potentials pi = de/dpi. At zero temperature they reduce to the Fermi 
energies Epi = yk Fi + M* 2 and the nonextensive statistical effects disappear. The meson fields are obtained as a 
solution of the field equations in mean field approximation and the related meson- nucleon couplings {g a , g u and g p ) 
are the free parameters of the model. In the following, they will be fixed to the parameters set marked as GM2 of 
Ref.fH. 

The thermodynamical quantities can be obtained from the thermodynamic potential in the standard way. More 
explicitly, the baryon pressure Pb and the energy density cb can be written as 

Pb = \T, j (S Xy [»?(*) +*?(*)] - \m 2 y - U(a) + \mtu 2 + \m 2 p p 2 , (18) 

i—n,p % 

** = 2 E / ^E*(k)[nJ(k)+nUk)} + \ m 2 a a 2 + U(a) + \m 2 ^ 2 + \ m 2 p p 2 . (19) 



It is important to observe that Eq.s ([T4|) . (TT8|) and (fl9|) apply to n\ = (rii) q rather than m itself, this is a direct 
consequence of the basic prescription related to the g-mcan expectation value in nonextensive statistics [20, [4y] (this 
recipe was not adopted in Ref.[42]). In addition, since all equations must be solved in a self-consistent way, the 
presence of nonextensive statistical effects influences the many-body interaction mediated by the meson fields. 

Especially in regime of low density and high temperature the contribution of the lightest mesons to the thermody- 
namical potential (and, consequently, to the other thermodynamical quantities) becomes relevant. As quoted before, 
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FIG. 1: The a meson field as a function of baryon chemical potential for different values of temperature (in units of MeV) and 



following Ref.s [54j, [55[ , we have included the contribution of pions considering them as a (quasi) ideal gas of nonex- 
tensive bosons with effective chemical potentials expressed in terms of the corresponding effective baryon chemical 
potentials. More explicitly, for tt + mesons we have p^+ — pc = Pp — Pn, where pc is the electric charge chemical 
potential. Thus, the corresponding effective pion chemical potential can be written as 



/4+ = Mp - /4 =tip~ Hn -2g P p, 



(20) 



where the last equivalence follows from Eq. (|17[) . Therefore, the p meson field couples to the total isospin density, 
which receives a contribution from nucleons and pions. 

Let us start our numerical investigation by considering the behavior of a, w and p meson fields at a fixed value 
Z/A = 0.4, for different values of temperature and nonextensive parameter q. Because meson fields have their source 
in the baryon and scalar density, which are very sensible to the behavior of the mean occupation number, all meson 
fields appear to be significantly changed in presence of nonextensive effects. 

In Fig. [TJ we show the a meson field as a function of the baryon chemical potential \Xb ■ It is interesting to observe 
that at lower pg, in presence of nonextensive effects, the value of the meson field is significantly increased for all 
values of temperature respect to the standard case, the other way round happens at higher ps- This important 
feature is due to the fact that, as already remarked in Section [Til for q > 1 and fixed baryon density (or fj,g), the 
(normalized) mean occupation function is enhanced at high values of its argument and depressed at low values. Being 
the argument of the mean occupation function Xi = /3(E? ~ /!*), in the integration over momentum (energy), at lower 
Pb (corresponding to lower values of the effective particle chemical potential p*) the enhanced Tsallis high energy 
tail weighs much more that at higher ps where depressed low energy effects prevail and the mean occupation number 
results to be bigger for the standard Fermi-Dirac statistics. Concerning the antiparticle contribution, the argument of 
Hi is Xi = f3(E* + p*) and the Tsallis enhancement at high energy tail is favored also at higher ps- At the same time, 
higher temperatures (where antiparticle contribution are more relevant) reduce the value of the argument of n% and 
Hi, favoring the extensive distribution. These effects are much more evident for the scalar density p$ (self-consistently 
related to the a meson field) where appears (rii) q and particle and antiparticle contributions are summed. The same 
effect involves also the nucleon effective mass M* — M — g a a, which becomes, respect to the standard case, smaller 
for lower values of pb and bigger for higher values, with very relevant consequences for the hadronic EOS 65]. 

In Fig. [2j we report the w meson field as a function of the baryon chemical potential for different values of 
temperature and q. In this case the situation is different from the a meson, because the uj field have its source in the 
baryon density ps where appears n, and particle and antiparticle contributions are subtracted. At lower temperatures 
(T = 60 MeV) , antiparticle contributions are negligible and we have a behavior similar (although less evident) to the 
a field. At higher temperatures (T = 120 MeV), the contributions of antiparticle increase and nonextensive effects 
vanish at higher ps- 

Finally, in Fig. [3] we report the behavior of the p meson field which depends from the isospin density (let us 
remember that we have fixed Z/A = 0.4). Similar arguments as done for the u> meson applies also in this case. The 
valuable increasing of its absolute value, also for weakly asymmetric nuclear matter, makes p meson very relevant in 
the hadronic EOS, especially at not too large ps- 
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FIG. 2: The u meson field as a function of the baryon chemical potential for different values of temperature and q. 
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FIG. 3: The p meson field as a function of baryon chemical potential for different values of temperature and q. 



In Fig. |U the total pressure P and energy density e are plotted as a function of fiB for different values of temperature 
and q. The different behavior from P and e reflects essentially the nonlinear combinations of the meson fields and the 
different functions under integration in Eq.s (|18[) and (|19j) . Concerning the pressure, we have that becomes stiffer by 
increasing the q parameter. On the other hand, the behavior of the energy density presents features very similar to 
the a field one. At low jUb, nonextensive effects make the energy density greater with respect to the standard case. At 
medium-high /i^, the standard (q = 1) component of the energy density becomes dominant, this effect is essentially 
due to the reduction of the a field for q > 1. The intersection point depends, naturally, on the physical parameters of 
the system. 



IV. NONEXTENSIVE QGP EQUATION OF STATE 



Concerning the nonextensive quark-gluon EOS, due to its simplicity, we adopt the MIT bag model [561 ]. In this 
model, quark matter is described as a gas of free quarks and all non-perturbative effects are simulated by the bag 
constant B which represents the pressure of the vacuum. 

Following this line, the pressure, energy density and baryon density for a relativistic Fermi gas of quarks in the 
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FIG. 4: Pressure (left panel) and energy density (right panel) versus baryon chemical potential for different values of temperature 
and q. 



framework of nonextensive statistics can be written, respectively, as 

d 3 k k 2 



^^ r°° c\ 3 k 

e q =7/ E I J^e f [nj(k)+nj(k)}+B, 



f=u,d' 



Pq 



2L 
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f=u,d 



d 3 k 

2^)3 

d 3 £; 

J2nf 



[n f (k)-n f (k)} 



(21) 
(22) 
(23) 



where the quark degeneracy for each flavor is 7/ = 6, ej = (k 2 + m 2 ?) 1 ' 2 , nf(k) and nf(k) are the g-deformed particle 
and antiparticle quark distributions 



n f (k) : 
n f (k) 



[l + (q-i)(ef{k)-»f)/T\i/ti-i) + l' 

1 

[l + (g-l)(e / (fc) + M/ )/r]i/( 9 -i) + l 



(24) 
(25) 



Similar expressions for the pressure and the energy density can be written for gluons treating them as a massless 
g-deformed Bose gas with zero chemical potential. Explicitly, we can calculate the nonextensive pressure P g and 
energy density e g for gluons as 



7g r d 3 fc k 

9 " 3 J (2tt) 3 [1 + (q - 1) k/T\9/(9-i) _ 1 
e g = 3 P g , 



(26) 

(27) 



with the gluon degeneracy factor j g = 16. In the limit q — > 1, one recovers the usual analytical expression: P g = 
8tt 2 /45T 4 . 

Let us note that, since one has to employ the fermion (boson) nonextensive distributions, the results are not 
analytical, even in the massless quark approximation. Hence a numerical evaluations of the integrals in Eq.s (|21l) - (|23p 
and (|26p must be performed. 

In Fig. [5j we report the total pressure as a function of the baryon chemical potential for massless quarks and 
gluons, for different values of q and at fixed value of Z/A = 0.4. The bag constant is set equal to B 1//4 =190 MeV. In 
presence of nonextensive effects, as in the case of hadronic phase, the pressure is significantly increased even for small 
deviations from standard statistics. 
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FIG. 5: Pressure of the quark-gluon phase as a function of baryon chemical potential for different values of temperature and 



V. THE HADRON TO QUARK-GLUON PHASE TRANSITION 



In this Section we are going to investigate the phase transition from hadronic matter to QGP at finite temperature 
and baryon chemical potential in the framework of nonextensive statistics. 

At this scope, we use the Gibbs formalism applied to systems where more than one conserved charge is present [57] . 
In fact, because we are going to describe the nuclear EOS, we have to require the global conservation of two "charges" : 
baryon number and electric charge. Each conserved charge has a conjugated chemical potential and the systems is 
described by two independent chemical potentials: ps and pc- The structure of the mixed phase is obtained by 
imposing the following Gibbs conditions for chemical potentials and pressure 



(H) _ (Q) 
p B — p B 



(H) _ (Q) 



P H (T,fi B ,^c)^P Q (T,fiB,fic) 



(28) 
(29) 



Therefore, at a given baryon density ps and at a given net electric charge density pc = Z/ApB, the chemical potentials 
Hb are p,c are univocally determined by the following equations 



Pb = (l-x)pf(r,MB,Mc)+XPi(r,MB,Mc), 

PC = (1 -X)Pc( T ^PB,^c)+XPc( T ^B,P-c), 



(30) 

(31) 



where p B and p c are, respectively, the net baryon and electric charge densities in the hadronic (H) and in the 
quark (Q) phase and \ is the fraction volume of quark-gluon matter in the mixed phase. In this way we can find out 
the phase coexistence region, for example, in the (T 7 pb) plane. We are particularly interested in the lower baryon 
density (baryon chemical potential) border, i.e. the first critical transition density p\ x (Mcr)i in order to check the 
possibility of reaching such conditions in a transient state during a heavy-ion collision at relativistic energies. 

In Fig. [HI we report the pressure as a function of baryon chemical potential; as before, we have set Z/A = 0.4 
and B 1 / 4 — 190 MeV. We can see that at T = 60 MeV, in the range of the considered ps, the system is in a 
pure hadronic phase even for the nonextensive index q = 1.1. At T = 90 and 120 MeV, we have both the first 
and the second transition critical density for the considered values of q. In presence of nonextensive effects, the 
values of the critical densities result to be sensibly reduced with respect to the standard case. This matter of fact 
is more evident in Fig. where we report the pressure at T = 90 MeV as a function of baryon density (in units of 
nuclear saturation density po = 0.153 fm~ 3 ) (left panel) and energy density (right panel). It is interesting to observe 
that pressure as a function of baryon density (or energy density) is stifier in the pure hadronic phase for q > 1 but 
appears a strong softening in the mixed phase. This feature results in significant changes in the incompressibility and 
may be particularly important in identifying the presence of nonextensive effects in high energy heavy ion collisions 
experiments. Related to this aspect, let us observe that possible indirect indications of a significative softening of the 
EOS at the energies reached at AGS have been discussed several times in the literature 39, 45, 58-60]. 

In Table U we report the critical baryon densities and baryon chemical potentials at the beginning (index /) and 
at the end of the mixed phase (index II) for different values of temperature and q. 
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FIG. 6: Pressure versus baryon chemical potential in the mixed phase for different values of temperature and q. 
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FIG. 7: Pressure as a function of baryon density (left panel) and energy density (right panel) in the mixed phase for different 
values of q. The temperature is fixed at T — 90 MeV. 



T = 60 MeV pj/po p'J/pp fjj [MeV] gg [MeV] 



q= 1.00 


5.75 


9.10 


1503 


1569 


q= 1.05 


5.56 


8.88 


1472 


1537 


q = 1.10 


5.33 


8.65 


1437 


1502 


T = 90 MeV 


pL/po 


pill Pa 


/4r [MeV] 


l& [MeV] 


g = 1.00 


3.41 


5.09 


1123 


1170 


g= 1.05 


2.77 


4.46 


1034 


1068 


g= 1.10 


1.91 


3.69 


916 


927 


T = 120 MeV 


pL/p0 


p'J/po 


/4 r [MeV] 


Mcr [MeV] 


g= 1.00 


0.45 


1.93 


588 


616 


q= 1.05 


0.20 


1.33 


383 


396 


g = 1.10 


0.08 


0.71 


184 


201 



TABLE I: Critical baryon densities and baryon chemical potentials at the beginning (index /) and at the end (index //) of the 
mixed phase for different values of temperature and nonextensive parameter q. 
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60 ' ' ' ' ' ' ' ' ' ' — ^^ ' ' ' ' ^^ Pb/Po 

2 4 6 8 

FIG. 8: Phase diagram T — ps for different values of q. The curves with index I and II indicate, respectively, the beginning 
and the end of the mixed phase. 



In Fig. El it is reported the phase diagram in the plane T — ps for different values of q. The curves labelled with the 
index / and // represent, respectively, the beginning and the end of the mixed phase. For q > 1, both the first and the 
second critical densities are sensibly reduced, even if the shape of the mixed phase is approximately the same. Related 
to this aspect, let us mention that the simplest version of the MIT bag model, considered in this investigation, appears 
to be not fully appropriate to describe a large range of temperature and density. To overcome this shortcoming, a 
phenomenological approach can therefore be based on a density or temperature dependent bag constant [45|, I61rro3j . 
Moreover, as discussed in the Introduction, in regime of high temperature and small baryon chemical potential the first 
order phase transition may end in a (second order) critical cndpoint with a smooth crossover. These features cannot 
be incorporated in the considered mean field approach. In our investigation, because we are focusing to nonextensive 
statistical effects on the nuclear EOS, instead of introducing additional parameterizations, we work with a fixed bag 
constant and limit our analysis to a restricted range of temperature and density, region of particular interest for high 
energy compressed nuclear matter experiments. 

Let us now explore in more details the variation of the first transition baryon density p cx as a function of different 
physical parameters. In Fig. [9] we report the dependence of p T cr as a function of Z/A for different values of q. It 
is interesting to note a significant reduction of p^ r in presence of nonextensive statistics; this effect is remarkable 
especially for T = 90 MeV and is essentially a consequence of the p meson field behavior in the hadronic phase. 

In Fig. [TU] (left panel), we show the first critical baryon density as a function of the bag constant for different 
values of nonextensive parameter q. Obviously, by increasing the bag constant we have a corresponding increasing of 
p^ r . However, this effect depends on the temperature and nonextensive parameter q. Finally, in the right panel of 
Fig. I10[ we show the variation of p^ T as a function of the nonextensive index q for different values of temperature and 
B X I A = 190 MeV. At T = 60 MeV, we can see only a little reduction in the first critical density also for large deviations 
from the standard statistics; on the other hand, the reduction becomes more pronounced at larger temperatures. 

VI. CONCLUSIONS 

To summarize, we have studied the main features of the nuclear EOS in the hadronic and quark-gluon phase 
and the possible formation of a consequent mixed phase in presence of nonextensive statistical effects. We have 
focused our investigation in regime of finite temperature and baryon chemical potential, reachable in high-energy 
heavy-ion collisions, for which the deconfinemcnt phase transition can be still considered of the first order. From a 
phenomenological point of view, the nonextensive index q is considered here as a free parameter, even if, actually 
should not be treated as such because, in principle, it should depend on the physical conditions generated in the 
reaction, on the fluctuation of the temperature and be related to microscopic quantities (such as, for example, the 
mean interparticle interaction length, the screening length and the collision frequency into the parton plasma). We 
have restricted our investigation for small deviations from the standard statistics and for values q > 1 because, as 
quoted in the Introduction, these values were obtained in several phenomenological studies in high energy heavy ion 
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FIG. 9: Variation of the first transition baryon density as a function of the net electric charge fraction Z/A and for different 
values of q. 
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FIG. 10: Variation of the first transition baryon density as a function of the bag constant (left panel) and nonextensive index 
q (right panel) for different temperatures. 



collisions. In this context, it is relevant to observe that by fitting experimental observable at q > 1, the temperature 
(or slop e) parameter T is usually minor of the one obtained in the standard Boltzmann-Gibbs statistics (q = 1) 
[24J |32|. This feature is also present in the considered nuclear equation of state because, at fixed energy per particle 
E/N, we obtain for q > 1 lower values of temperature respect to the standard case. Moreover, let us remember that, 
in the diffusional approximation, a value q > 1 implies the presence of a superdiffusion among the constituent particles 
(the mean square displacement obeys to a power law behavior (x 2 ) oc t a , with a > 1) [64J . 

In the first part of the work, we have investigated the hadronic equation of state and the role played by the 
meson fields in the framework of a relativistic mean field model which contains the basic prescriptions of nonextensive 
statistical mechanics. We have shown that, also in presence of small deviations from standard Boltzmann-Gibbs 
statistics, the meson fields and, consequently, the EOS appear to be sensibly modified. In the second part, we have 
analyzed the QGP proprieties using the MIT Bag model and also in this case the EOS becomes stiffer in presence of 
nonextensive effects. Finally, we have studied the proprieties of the phase transition from hadronic matter to QGP 
and the formation of a relative mixed phase by requiring the Gibbs conditions on the global conservation of baryon 
number and electric charge fraction. We have seen that nonextensive effects play a crucial role in the deconfinement 
phase transition. Moreover, although pressure as a function of baryon density is stiffer in the hadronic phase, we have 
shown that a strong softening in the mixed phase takes place in presence of nonextensive statistics. Such a behavior 
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implies an abruptly variation in the incompressibility and could be considered as a signal of nonextensive statistical 
effects in high energy heavy ion collisions. 
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